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Abstract. We analyze the quantum regime of the dynamical backaction cooling of a 
mechanical resonator assisted by a driven harmonic oscillator (cavity) . Our treatment 

applies to both optomechanical and electromechanical realizations and includes the 
effect of thermal noise in the driven oscillator. In the perturbative case, we derive 
the corresponding motional master equation using the Nakajima-Zwanzig formalism 
and calculate the corresponding output spcc;trum for the optomechanical case. Then 
we analyze the strong optomechanical coupling regime in the limit of small cavity 
linewidth. Finally we consider the steady state covariance matrix of the two coupled 
oscillators for arbitrary input power and obtain an analytical expression for the final 
mechanical occupancy. This is used to optimize the drive's detuning and input power 
for an experimentally relevant range of parameters that includes the "ground state 
cooling" regime. 
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1. Introduction 

Recent progress in the emerging fields of nanoelectromeclianical [H |2| and 
optomeclianical systems [31 H] promises to enable quantum limited control of a single 
macroscopic mechanical degree of freedom [Sj IH]. This is relevant in the context of high 
precision measurements [Tj-[in] like single spin magnetic resonance force microscopy 
[TT| [T2] (MRFM) and for fundamental studies of the quantum to classical transition 
[T3l [HI [15] . A paradigmatic goal which has triggered a surge of activity is to prepare the 
eingenmode associated to an ultra long lived mechanical resonance (angular frequency 

and Q- value Qm) in its quantum ground state with high fidelity [IS]-[2S]- The 
considerable difficulty to achieve the desired combination ksT -C hum, Qm ^ 1 with 
state of the art micro-fabrication and cryogenic techniques [2] has naturally motivated 
ideas to use cooling schemes analogous to the laser-cooling of atoms [29]-[32]. 

In these schemes the mechanical resonator's displacement is coupled parametrically 
to an auxiliary high frequency bosonic or fermionic resonator (pseudospin) that can act 
as a "cooler" [33]. To drive the latter while monitoring its output allows detection of 
the mechanical displacement. Naturally there is a back-action force associated to this 
measurement process [311 EHl [36]. Due to the dissipative dynamics of the cooler for 
a negative detuning of the drive this force becomes anti-correlated with the Brownian 
motion resulting in net cooling. In turn the quantum fluctuations of the cooler — which 
in the atomic laser cooling manifest in the inherent stochastic nature of the spontaneous 
photon emissions — result in a quantum noise spectrum for this backaction force that 
sets a fundamental lower bound for the final temperature [301 EI] • Thus the structured 
reservoir afforded by the driven cooler and its environment provides an effective thermal 
bath for the mechanical resonator. The concomitant absorptions of motional quanta 
(cooling) correspond to Raman scattering processes in which a drive quanta is up 
converted, while emission events (heating) are associated to Raman processes in which 
a drive quanta is down converted [cf. Fig. [l|^b)]. 

A host of concrete realizations of the above generic scenario have been discussed 
in the literature. These range from electronic or electrical devices where the cooler 
is provided by a (superconducting) single electron transistor [TH], a Cooper-pair box 
[37] . an LC-circuit [3j, or a quantum dot to optomechanical systems where 

this role is played by an optical cavity mode [1]. In the latter systems, which are 
equivalent to a Fabry-Perot with a moving mirror [cf. Fig. [l|^a)], the optical field 
couples parametrically to the mechanical motion via radiation pressure. This effect 
has already been thoroughly demonstrated experimentally [39] and harnessed to provide 
appreciable cooling, albeit in the classical regime [201 EH ESI [261 EZ] ■ Recent experiments 
involve both optical and microwave cavities. In turn completely analogous effects 
have been shown using a capacitively coupled radio frequency LC-circuit [21] and a 
superconducting microwave cavity [28j. These systems share the common feature that 
the cooler is well approximated by a single harmonic oscillator which henceforth will 
be referred to as the "cavity". While for optical cavities the vacuum constitutes an 
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Figure 1. (a) Fabry-Perot equivalent of a mechanical eigenmode (frequency cjm/27r 
and Q- value Qm) coupled to an optical cavity mode, (b) Level diagram of the 
two modes in a "shifted" representation for perturbative optomechanical coupling 
T]. Raman scattering processes induced by the latter can decrease (solid lines) or 
increase (dashed) the mechanical eigenmode's quantum number n (a is the steady 
state amplitude in the cavity mode and |0)j,, |l)p,... its Fock states). [Reproduced 
from I. Wilson-Rae, N. Nooshi, W. Zwerger, and T.J. Kippenberg. Theory of ground 
state cooling of a mechanical oscillator using dynamical backaction. Phys. Rev. Lett., 
99:093901, 2007. "Copyright (2007) by the American Physical Society."] 

excellent approximation for the input when the drive is switched off, in the case of radio 
and microwave frequencies thermal noise in the cavity needs to be taken into account. 
A quantum treatment of the corresponding temperature limits has already been given 
which predicts that ground state cooling is possible when the mechanical oscillation 
frequency is larger than the cavity's linewidth k |101 HU, Il2] . Here we provide a rigorous 
analysis of the cooling dynamics that drives the mechanical resonator mode to a thermal 
state with a well defined effective final temperature that for a finite Qm is imprinted 
on the cavity's output. This is done in Sec. [3] where we obtain the motional master 
equation and the corresponding output spectrum. 

The above picture is only valid for perturbative optomechanical coupling and 
thus breaks down when the resulting cooling rate becomes comparable to the cavity 
linewidth k, or to the mechanical oscillation frequency [H] Um- In the Doppler limit 
ujm <^ the system then approaches the parametric instability or settles into a regime 
where backaction effects are purely diffusive with no net cooling. On the other hand 
in the resolved sideband limit k <^ Um as the optomechanical coupling exceeds the 
cavity linewidth the system enters into a strong coupling regime in which the motion 
hybridizes with the cavity fluctuations. The ensuing optomechanical normal modes are 
then cooled simultaneously. This phenomenon is analyzed in Sec. |4]in the limit of small 
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cavity linewidth where we find that the dynamics of each of the normal modes can be 
described by a master equation analogous to the one valid in the perturbative regime 
with a cooling rate given by half the cavity linewidth. Finally, in the same Section we 
derive an analytical expression for the final (steady state) average mechanical occupancy 
(phonon number) valid for arbitrary optomechanical coupling and use it to optimize the 
parameters of the drive. 

2. Optomechanical Master Equation 

In optomechanical systems (and their electromechanical analogues) the cavity resonant 
frequency Up depends inversely on a characteristic length that is modified by the 
mechanical resonator's displacement. The fact that Up 2> allows for an adiabatic 
treatment of this effect which results in the aforementioned parametric coupling. In 
general the leading contribution to the latter is linear in the mechanical displacement 
— but situations in which it is instead quadratic can be readily engineered ^21 [26] . The 
desired cooling dynamics is induced by a slightly detuned electromagnetic drive (angular 
frequency ujl), which for optomechanical systems corresponds to an incident laser and 
in the electromechanical case is afforded by a suitable external AC voltage. Thus the 
Hamiltonian describing the coupled system (in a rotating frame at ul) is given by |13] 

H' = -^A'^ajop + hr]ujmalap{am + aj„) + h^{ap + aj) + huj^al^a„^ . (1) 

Here (am) is the annihilation operator for the electromagnetic (mechanical) oscillator 
and A'^ is the detuning of the drive from Up. Here we have defined f2 = 2^ PK^y^jfiuJi^^ 
where P is the input power of the drive and is the cavity decay rate into the associated 
outgoing electromagnetic modes. The dimensionless parameter 77, characterizing the 
non-linear coupling between the cavity and the mechanical resonator, is given by 
ri = {ujp/uJm){xo/ L) where xq is the zero point motion of the mechanical resonator 
mode and the characteristic length L depends on the physical realization. In the 
optomechanical case it corresponds to an effective optical cavity length while for 
electromechanical realizations [21] L = 2dCtot/Cc, where Cc oc 1/d is the dynamical 
capacitance, d is the distance between the corresponding electrodes, and Ctot is the 
total capacitance. 

We treat the losses induced by the electromagnetic and mechanical baths within 
the rotating wave Born-Markov approximation using the standard Lindblad form 
Liouvillians [H]. It is important to note that the validity of a rotating wave 
approximation (RWA) in the environmental coupling responsible for the mechanical 
losses only amounts to Qm ^ 1 provided the optomechanical coupling is weak enough 
that there is no appreciable mixing between the annihilation and creation operators 
of the modes (cf. Sec. Clearly if the latter is not satisfied the usual RWA will 
result in the unwarranted neglection of resonant terms. This can be borne out 
by comparing the corresponding displacement spectra and results in the condition 
r7|a;|a;TO -C max{y/uJ^, Um}- Henceforth we focus on parameters that satisfy it which 
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include the most relevant regimes for cooling and ensure that the system doe not 
approach the instability for any detuning. Thus the evolution for the density matrix of 
the resonator-cavity system reads 

+ ^[n{ujp) + l]{2appal - alapp - palap) 
+ ^n{uJm){'2a'l,pam - a„,al^p - pa^a}^^ 

+ + l](2«mPaL - «m«mP " pamOm) • (2) 

The total cavity decay rate k has two contributions: (i) the rate at which photons are 
lost from the open port (where the driving field comes in) given by and (ii) the 
"internal loss" rate k — Kex due to the other losses of the electromagnetic resonator 
(i.e absorption inside the dielectric, scattering into other modes, etc.). Naturally 
the thermal noise is determined by the Bose number niiS). At room temperature 
and for optical frequencies n(u;p) is negligible, however for the much lower radio and 
microwave frequencies characterizing electromechanical setups this quantity can be 
comparable to the final mechanical occupancies achieved. Similarly 7^ = oJmlQm is 
the mechanical resonator's natural linewidth and n{(jJm) its mean occupation number at 
thermal equilibrium (i.e. in the absence of the drive). 

To study the cooling process it proves useful to apply a canonical transformation of 
the form — Op + a, — a™ + /? with the amplitudes a, [3 chosen so that the linear 
terms in the transformed Liouvillian cancel out. This condition leads to the following 
coupled equations for the amplitudes 

- (A'^ + i'^)a + 7]u:ma{l3 + /?*) = , 
(cUm - i^)/? + r^t^mlap = . (3) 

We assume 77 <^ 1 and that the mechanical dissipation rate 7^ is much smaller than 
uOm- To lowest order in the small parameters r] and 1/Qm we obtain a ~ fi/(2A'j^ + ik) 
and (3 ~ — ?7|ap. Here |ap is the steady state occupancy of the cavity in the absence 
of optomechanical coupling and [3 is the static shift of the mechanical amplitude due to 
the radiation pressure. The normal coordinates after the transformation are shifted so 
that the new amplitudes correspond to the deviation from the steady state equilibrium 
position. This transformation leaves the dissipative part of the Liouvillian invariant and 
transforms the Hamiltonian into: 

H = —hAio^pap + fkUmal^am + h-qurnia^cip + tt*«p + aaj,)(am + aL) ■> (4) 

where we have introduced the effective detuning A'^ + 2rf\a\'^uJm A^. It is interesting 
to note that if the bosonic cooler is replaced by a fermionic one so that — > cr_ we 
obtain a Hamiltonian that resembles the one describing a trapped ion in the Lamb-Dicke 
regime. As a result for perturbative rj (as analyzed in the next Section) the cooling cycle 
[cf. Fig. W(h)] becomes analogous to the Lamb-Dicke regime of atomic laser-cooling [lO] . 
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3. Perturbative Cooling 

3.1. Master equation for mechanical motion 

We first focus on tlie regime in wfiicfi tlie input laser power P is low enough so that the 
time scales over which the populations of the mechanical resonator's Fock states evolve 
(leading to cooling or heating) are much slower than those associated to the losses of the 
cavity and to the free mechanical frequency. As will become clear below this requires 

Here we also assume [n{ujm) + l]7m ^ i^, '^m which must hold to allow for 
appreciable cooling, and <^ \ (in optomechanical realizations ?7 < 10^^). Hence 
the electromagnetic degrees of freedom can be treated as a structured environment 
that affects the mechanical motion perturbatively. Along these lines the latter can be 
described by a Markovian master equation for its reduced density matrix |44j . 

To derive it we take the optomechanical master equation in the shifted 
representation, transform to an interaction picture for the resonator mode and 
adiabatically eliminate the cavity using the Nakajima-Zwanzig formalism [l5l HU B6] . 
The optomechanical coupling and the mechanical losses are treated perturbatively. We 
define the projection 

Pp = Tr4p}®pf) Q = (5) 

with 



Pf^ = - 



oo 

^5: 



Up 



^_np + 1 



n){n\p, (6) 



where Up = n{ujp) is the Bose number, and introduce the formal parameter ( such that 

£(t) = C% + C/:i(C't) + A, (7) 

with 

Cop = i[ALaJap, p] + -n{ujp)[2alpap - a^ajp - pflpaj] 

+ -[n{ujp) + l][2appaj - aja^p - pajop] , (8) 

Ci (C't) = e'"'"^'*/:S+^ + e-'"'"^'*£S-) , (9) 

C[^^p = -i [i]uJm [flpflp + a*ap + aal] a]^, p] , (10) 

= Vl^m [alap + + aal] ^rn^ p\ ' (^^) 
^2 = ^n{uJm){2al^pa^ - a^a^^p - pa^al^) 

+ ^[n{uj^) + l]{2a^pal^ - al^a^p - pal^aj . (12) 

We note that Vp is a stationary state of Cq for any p implying that CqV = while 
VCq = follows from trace preservation, so that we have 

QCoQ = Co, VCoV = QCoV = VCoQ = 0. (13) 
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As will emerge from our derivation the basic idea is that the rates for cooling and heating 
set the relevant time scale (zeroth order in l/Q which is widely separated from the 
mechanical oscillation period 2tt /ujm and the cavity ring down (order 1/C^)- In fact, 
the asymptotic expansion for ^ ^ oo pursued below amounts to a controlled expansion 
in the ratio between the fast and the slow timescales. Given that we are interested in 
the behaviour as t — > cxd, the initial condition is immaterial and we choose for simplicity 
one in the P- manifold so that Qpo = 0. Subsequently by explicitly integrating the 
differential equation for the irrelevant part Qp, we obtain a closed equation for the 
relevant part Vp: 



Vf> = VC{t)Vp + VC{t) / drT^ 



Jo'dr'Q£(r')S 



^dr" QC(t")Q 



C{T)Vp{r) , (14) 



where 7^ (71) is the time-ordering (anti-time-ordering) operator. For the purpose of 
analyzing the asymptotic limit ( 

'1 



J^dr'QC{T')C 



■£dT"QCir")C 



OO of Eq. (|14|) we have 

O 



,CQCoQ{t~T) 



1 



(15) 



which can be understood considering the corresponding Laplace transforms. 
Substitution of Eqs. (13), (15) into Eq. (14) and the change of variables r' = C^(t — r) 
then yields 



(16) 



Vp = P[CC,iC^t) + C2] Vp 

+ VC^{Ch)Q / dr'e^°^'Q/:i(C'i-r)Pp(t-r7C') + ... 
Jo 

Here we have also used = Q, = V, and VC2 = C2V. The leading order in 1/C 
of the evolution over the aforementioned relevant time scale will be determined by the 
limit as C — ^ CX3 of the Laplace transform of the above. In the time domain all the fast 
rotating terms (frequencies of order C'^uJm) drop out and Eq. (16) reduces to 



Vp = VC2VP + 



poo 

VC[^^Q / dr'e(*"™+^")"'Q/:l"^Pp + H.c. 
Jo 



We note that as V 
^op = 



'0 

= lim^^oo e^"* we have 
^ Qp = 0. 



(17) 



:i^ 



It follows that the restriction of Cq to the Q-manifold has only eigenvalues with negative 
real parts which allows us to establish 



POO 

Q / dr'e(*-'-+^°)^'Q= Q(-^a;„-/:o)~'Q• 
Jo 



(19) 



We focus on the parameter regime where |ap ^ Up. The behavior of the cavity 
correlations that determine the second term in Eq. (17) implies that in this regime 
contributions arising from the cubic term in C^^^ are negligible compared to those 
generated by the quadratic term — for Hp = the contributions of the former are 
higher order in 7]. The conditions that warrant this linearization of C^^^ for the case 
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will be discussed further in the next Subsection. The restriction of the quadratic 



term to the P-manifold vanishes and we obtain 
Trp !^VC[^^ Q dr' e(^'^'"+^«)"' QC['^Vp^ 



^ - Y {G{^rn, rij,) [al^, a^fj] - G'*(-^™, n.p) [a|„, /xa J } . (20) 
Here we have used Eqs. (j5]), (10), (11), and introduced the optomechanical coupling 

gm = '2r]\a\um , (21) 

the reduced density matrix for the mechanical mode /i = Trp{p}, and the cavity- 
quadratures' correlations 



G{Lj,np) 



dre*-^Trp{Xp(0)e^«^Xp(0)pf)} , 



(22) 



with Xp{0) = {a* Up + aal)/V2\a\. If we substitute Eqs. (|5|, ([l2|, (|20|) into (|l7|, trace 
out the cavity, and rearrange we finally obtain the following master equation 

/i = -i[{uj.m + A.^)al^am, l4 

+ ^{lm[n{ujm) + 1] + v4_(np)}(2amyua]„ - al^amj^ - na\^am) 



(23) 



where the cooling (heating) rate A_{np) [A^{np)] and the mechanical frequency shift 
Am induced by the optomechanical coupling are given by 

^t(^p) = 9rn^{Gi±UJm, rip)} , 



Am = ^Q{G{uJm, np) - G{-UJm, ^p)} • 



(24) 
(25) 



It is straightforward to calculate the necessary two-time correlations using the quantum 
regression theorem given the steady state moments: 

(a^aj) = rip + 1 , (aja^) = Up , (apOp) = , (26) 
and that the evolution of the mean amplitude reads 

{hp) = {zAl - k/2) {ap) . (27) 
Thus we obtain 

G{uj, Up) = [G{uj, 0) + G* (to, 0)] Up + G{uj, 0) , (28) 

with 

G{uj,0)= _ \ ^ , ^ (29) 



-2i(uj + AL) + K 



which substituted into Eqs. (|24j), (25) leads to 

A^irip) = [A_{0) + A+{0)]np + A^{0) 



A 



At. - UJ„ 



+ 



Al + UJr, 



4(Ai -U:mY + 4( Ai ^UJmY^^^^\ 



(30) 
(31) 
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where the corresponding rates in the absence of thermal noise in the driven cavity are 
given by: 



The above can be related to the input power and the frequency of the drive via: 



^- = 2^^™^/n.,(Ai + .V4)' ^^^^ 



The master Eq. (23) generalizes the one obtained in Ref. [ID] by including thermal 
noise in the cavity input. As shown below, this effect is significant for determining 
the ultimate limit to which the mechanical resonator can be cooled for ratios Up/um 
like the ones that characterize electromechanical realizations. Equation (23) has 
as its steady state a thermal state which defines the effective final temperature to 
which the mechanical resonator is cooled. The corresponding final occupancy reads 
[r = A_(0)-A+(0)]: 

where ^/ is the quantum backaction contribution derived in Refs. JlHlllI], namely 

~ ^ (Al + UJmf .gg. 

^ 4Alu;„ 

We note that in the "unshifted" representation there is in addition a coherent shift of 
the resonator's normal coordinate so that 



{aln^J =nf + T] 



Pk, 



(37) 



If we now consider the appreciable cooling limit F ^ 7^ and minimize with respect to 
the detuning we obtain 



min{n^} ^ ^ ^^^^ ""''^ + + (n^ + [ ^1 1 + ^ ~ 1 ] , (3J 

for the optimal detuning A^^* = ^^/uJ^^^^^~^^^/4:. The first term corresponds to the 
"linear cooling" limited by thermal noise. In turn, the second term shows that the final 
occupancy is necessarily bounded by the equilibrium thermal occupancy of the cooler. 
Finally, the last term for Up = corresponds to the fundamental temperature limit 
imposed by the quantum backaction which in the Doppler regime k <^ Um reduces to 
n/iujrn, precluding ground state cooling, while in the RSB regime it yields K^/16a;^ 
corresponding to occupancies well below unity jlQl HI]. We note that for a given 
K/iOm < v^32 occupancies below unity are only attainable within a finite detuning 
window \Al + 3a;m| < ^/8ui^^^-^~K^j4:, as illustrated in Figure [2] 
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3.2. Output power spectrum and temperature measurement 

The spectrum of the cavity output constitutes a crucial observable to understand the 
backaction cooling. It allows to visualize the cooling cycle as a frequency up-conversion 
process and for optomechanical realizations it provides an efficient way to measure 
the final temperature ^0]. We focus on the latter for which Up = and consider 
the experimentally relevant case in which the output is measured in the same modes 
in which the coherent laser drive is fed [22] • To calculate its spectrum we apply the 
standard input-output formalism [H] and treat the parameters ri\a\, rj perturbatively 
along the lines of the previous Subsection. For this purpose it proves useful to consider 
the quantum Langevin equation for Up associated to the Liouvillian ([T]), namely 

dp = (^iAl - flp - ir]ujm{a^ + + a) + y/i^Sam{t) + V - /«cx&in(i) , (39) 

where 6ai^ and b^^ correspond to the vacuum noise associated, respectively, to the laser 
mode and to the other cavity losses. Here we use the shifted representation [cf. (|4])] and 
also take into account that Hamiltonian ([T| presupposes the standard time-dependent 
canonical transformation that maps the coherent input state associated to the laser into 
a classical field so that ain(t) —>■ 5ciin(^) + {(iin{t))- If we assume that the term linear in r) 
is a specified function of time and that the solution a^\t) for rj = is known, Eq. (39) 
can be formally integrated to obtain 

a^t) = - ^V^m fdr e('^-''/2)(*-) [aUr) + ajr)] [ap(r) + a] . (40) 





This integral equation can be iterated to generate the following Dyson series type result 



apit) = af (t) + f^{-t^uj^)- f dr„ Tdr^^i . . . H dr, e^^^^^^/^^^'-^^^ [a|„(r„) 

Jo Jo Jo 

+a^(r„)] . . . [alin) + a^{n)] [a^^\n) + a] . (41) 
In turn, the output spectrum is given by 



1 r+oo 

S{u;) = - dr e-(----)^(aL(t + r)a,uMss , (42) 
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which in the shifted representation reads 
1 

X 



S{u;) 



dr e 



(t + r) + (5aL(t + r) + v^a* + (ain(O))^ 



/K;;;ap(t) + 5ajj^(t) + ./k;;;;^ + (ai^(o)) )ss- (43) 

We now seek the lowest non-trivial order in 77. The output flout has a c number part 
arising from the classical drive and the cavity steady state amplitude and an operator 
part corresponding to the fluctuations. Equation (43) directly implies that terms 



involving the c numbers only contribute to the "main line" oc 5{uj — ujl). Furthermore, 
as in the shifted representation the steady state of the electromagnetic modes is the 



vacuum |0) it follows from Eqs. (41), (43) that the contribution of the cross term with 
the operator part — which vanishes if the cubic term in Eq. (|4]) is omitted — is at 
most higher order in rj and can be neglected relative to the contribution bilinear in 
the c numbers. The latter corresponds to the classical reflection coefficient that is 
straightforward to obtain considering the near-resonant scattering into modes (real or 
fictitious) responsible for the other losses k — k^^. Thus we arrive at 



S{u;) 



P 

hjjT 



1 - 



X 



Ai + /tV4 

+00 

dr e"^('^~^^)^ 



5{ij - lul) + 



in 



t + T 



dr{ e*- 



aL rJ + a. 



r. 



X 



dri e 



{iAL-K/2)it-Tj) r t 



[alin) + a„(ri)] 



(44) 



ss 



where we have used a^\t)\0) 



5ain|0) = and the correction is higher order in rj 



for all frequencies. To calculate the steady state two time average in Eq. (44) we 



adopt (as in the previous Subsection) an interaction picture for the resonator mode 



[i.e. am{t) 



-iulmt 



amit)] and make the substitutions t + r — r{ — r(, t + ri — > ri in 



the time integrals. In this representation the mechanical mode operators evolve slowly 
compared with 1/k so that in line with our derivation of a Markovian master equation 
for the mechanical motion we can factor them out of the time integrals whose upper 
limit can be extended to infinity. This Markovian approximation yields 

X [alin) + a„(ri)] \ ^ 2( [G*(^^, 0)e^--(*+-)aJ„(t + r) + G*{-u^, 0) 
/ ss 

xe-^--(*+-)a^(t + r)] [G(-cu^, 0)e^-™*aj„(t) + G(a;^, 0)e-^--*a„(t)] )m 



where we have used Eq. (28) with Up = 0. The two time averages of the mechanical 



mode operators can now be calculated from the master equation (23) using the quantum 
regression theorem. The needed one time averages satisfy 

7m + 



{am) 



iAm + 



{am) 



(«Lom)ss = nj , (a^aDss = ^/ + 1 , (a^«m)ss = {a^al^) 



ss 



0. 



(46) 
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Finally, a straightforward calculation leads from Eqs. (44), (45), (46), and (29) to the 
spectrum already given in Ref. 

iA^Uf 



P 



A| + /t74 
{uf + 1) 



EO], namely 

6{uj 



Toff 
2 



KTC 



Tcff 
2 



(47) 



7m + r which is the total dissipation rate for the 



where we have introduced 7efr 
mechanical mode in the presence of the drive that determines the linewidths of the 
motional sidebands peaked at ujl ±ujm- One should note that the above corresponds to 
photon rate per unit frequency and that the weight of the motional sidebands relative 
to the main line is of order rf'{2nf + 1) [77^ < 10~® for typical systems]. Here unlike the 
case of atomic laser cooling [301 EI] the mechanical dissipation induces an asymmetry 
in the weights of the sidebands that allows to retrieve the final temperature directly 
from the steady state. The "blue" sideband weighted by = ^A^Uf corresponds to 
the up-converted photons (anti-Stokes scattering) responsible for the cooling while the 
"red" sideband weighted by = + 1) corresponds to the down-converted 

photons (Stokes scattering) that result in heating. 

It is interesting to note that the formalism of this Section does not presuppose 
linearizing around the steady state and neglecting the cubic term in Hamiltonian 
^ accordingly, but rather the validity of such treatment emerges from a controlled 
procedure that would allow to incorporate the necessary corrections if the intrinsic non- 
linearity rj were larger. A straightforward self-consistency criterion is to compare the 
steady state fluctuation of the cubic term with that of the quadratic one. An heuristic 
estimate of their ratio can be extracted from the total weight of the motional sidebands 
in Eq. (47) which together with the analysis in Subsec. 3.1 implies that the cubic term 
can be neglected provided the conditions rf{2nf + l)c<j^ ^ and |ap ^ rf are satisfied 
— for the typical parameters in cooling experiments these are always met. 



4. Linearized theory for coupled optomechanical modes 

4-1. Small cavity linewidth limit n Qm 

The treatment in the previous section is only applicable when the cooling rate A^ given 



by Eqs. (30), (31), (32) is much smaller than k (note that the heating rate A^ is always 
bounded by A_ for negative detuning). This condition for arbitrary negative detuning 
results in g'^ <^ so that it follows (as expected) that the motional master equation is 
only warranted for small enough optomechanical coupling. In the Doppler regime when 
this is violated the aforementioned condition g^ <^ maxj^/cJ^, Um} underpinning the 
RWA for the mechanical losses will also fail. In contrast in the resolved sideband regime 
(RSB) relevant for ground state cooling n -C implies that there is a wide parameter 
range of interest in which Eq. ^ remains valid while Eq. (23 ) fails. Here we consider this 
RSB regime beyond perturbative optomechanical coupling. As gm becomes comparable 
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to K it becomes necessary to follow the coupled dynamics of both modes as described 
by Eq. ^ which for > k/2 exhibits normal mode splitting. Though for a Gaussian 
initial condition the approach to the steady state is always amenable to a straightforward 
description, in the intermediate regime gm ^ there will be no simple analog of Eq. (23) 
that allows to visualize the cooling process in terms of phonon jumps. In turn, deep in 
the strong coupling regime but away from the instability, i.e. for n <^ gm ^ ^m, the 
dynamics can be described by two decoupled master equations for the optomechanical 
normal modes analogous to Eq. (23). 

Within the latter parameter range it is permissible to start from a Hamiltonian 
description including the optomechanical coupling and treat the losses (as described by 
7m, perturbatively. We focus on the "resonant case" — = ojm — which in the RSB 
regime can be shown to be optimal for minimizing the final occupancy — and consider 
the canonical transformation that diagonalizes Hamiltonian (|4]) for 77 = with (^m 7^ 
[after performing the convenient rotation Op {a/\a\)aj^. The latter is given by 

1 




'"m/p 



where the eigenfrequencies of the normal modes read 



(4J 



(49) 



If we now consider the expansion of the above in the small parameter gm/^m the 
zeroth order of Eq. (49) yields a splitting given by gm while Eq. (48) reduces to the 
transformation that diagonalizes the rotating wave part of Hamiltonian (|4]) — that 
results from neglecting the terms that involve a^ap, a\^a]p. The latter transformation 
does not mix the annihilation and creation operators. For small but finite gm/^m there 
will be small admixtures that for the purpose of analyzing the cooling dynamics will 
only be relevant insofar as they give rise to qualitatively new terms in the dissipative 
part of the Liouvillian — otherwise they can be shown to result in contributions of 
relative order {gm/^mY for all values of the other parameters. Hence we have 



^m/p 



1 I 1 9m t I 9m 

—j=a+ ± —j=a^ -= — a\ ± — -= — 

y/2 y/2 Ay/2uJm Ay/2uJr, 



(50) 



To proceed we: (i) apply the transformation given by Eq. (48) to the master 
equation ([2]), (ii) transform the result to an interaction picture with respect to the 
(now diagonal) Hamiltonian, (iii) neglect all the resulting fast rotating terms which are 
cx e^^*'^""* or oc e^*^™* (up to corrections higher order in gm/^m) and (iv) expand to 
lowest order in gm/^m following the aforementioned "qualitative" criterion. Naturally 
(iii) relies on the small cavity linewidth condition n <^ g^. Thus we obtain 



P 



E 



7r7 



+ 



+ -n[uj. 



pj 



1^9'L 

m 



2a 



In 



1] I ^2agpa^ — a^a^p — pO'l^i^ ■ 



(51) 
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The only corrections in {gm/^rnf' appear in the first term and correspond to the "high 
power hmit" of the heating induced by the quantum backaction of the cavity. Analogous 
heating terms oc 7™ are neglected given that they are comparable to corrections to the 



RWA treatment of the mechanical dissipation. It follows from Eq. (51) that the losses 



do not couple the normal modes (annihilation operators a±) so that the steady state is 
given by the tensor product of thermal states for each of them characterized by average 
occupancies 

2 



aj_a_|_) 



ss 



(52) 



7m + 

(where we neglect higher order terms in gm/i^) to which the average phonon number 



converges with a cooling rate {k + 7m)/2. Hence Eqs. (50), (52) finally yield 



7„n(u;„) + Kn{ujp) + 



7m + K 

K 



+ rip + 



9' 
16u;2 

m 



(53) 



where the corrections are higher order in the small parameters {gm/^mY) Im/ 1^- The 
first term corresponds to the heating associated to the mechanical dissipation and can 



be identified with the corresponding term in Eq. (35) showing a saturation of the usual 



linear cooling law. Similarly, the second and third terms can be identified with the 



corresponding contributions in Eq. (35) arising from the thermal noise in the cavity and 



the quantum backaction. Thus ground state cooling requires ksT /hQm ^ k ^ ujm and 
Up < ksT. We note that Eqs. (50), (52) imply that (a^)ss = so that the reduced 



state for the motion is also thermal in the small k limit. 



If we now compare Eqs. (38) and (53) and consider minimizing them with respect 



to gm within their respective ranges of validity, heuristic considerations imply that the 
following formula is expected to always constitute a lower bound for the final mechanical 
occupancy optimized with respect to the parameters of the drive 



7mn{UJm) 



'1 + 



1 



(54) 



4-2. Final occupancy for arbitrary ratio gm/ 1^ and optimal parameters 



The approximate expressions (35), (53) for the final mechanical occupancy that we 
have derived in the limits g-m ^ 1^ and gm ^ 1^ provide a basic understanding of the 
requirements for ground state cooling and the expected order of magnitude for the 
optimum. Notwithstanding they have the drawback that they fail to settle which is 
the optimal input power as minimization with respect to gm shifts this variable away 
from the domain where they are valid. In addition given the experimental progress 
towards achieving ultra-cold states in these systems it is clear that precise quantitative 
predictions of the steady state for a given input are highly desirable. To this effect 
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we complement the above analysis by deriving an analytical expression for Uf valid for 
arbitrary values of the ratio gm/f^- The optomechanical master equation (|2| directly 
implies that when the cubic non-linearity in Hamiltonian ^ is neglected the time 
evolution for the ten independent second order moments that determine the covariance 
matrix is given by a linear system of ODEs. This is determined by 







(55) 






(56) 


d / t \ _ 




(57) 






(58) 


d / 2 \ 




(59) 



(60) 



dt 

and their Hermitian conjugates. The exact solution of the linear system of equations 
that results for the steady state covariance matrix then yields an analytical formula 
for the final occupancy of the mechanical resonator that is a sum of three independent 
contributions 



n 



Here n 



(m) 



/ 

,{P) 



n 



(m) 
f 



(0) 



(61) 



J and arise, respectively, from the mechanical dissipation and the thermal 
noise in the cavity input, and n^p corresponds to the heating induced by the quantum 
backaction exerted by the cavity. The expressions for these different contributions are 
rather unwieldy for arbitrary mechanical dissipation 7^, but it is simple to realize that 
for appreciable cooling n/ -C n{u!m) to be possible the conditions 7^ <^ k, gm, '^m are 
needed. Hence in this regime of interest a good approximation for is obtained if one 
takes 7m keeping ^.mniio^ finite. In this limit we obtain 



n 



(0) 

/ 

(m) 



+ A.)V- 



+ 



Ti \ 

f KAig^ul \ 4 



K 



{Um + Al)' + Y 



2 



U). 



A \ 2 ^ 

Al + -T 



K 

~2 



K 



n 



(p) 



Up R 



Al {u^Al + uj^Al + glAl/2 + ulgl) + {2u;mAl + ^^Ai/2 
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(b) (1.20 - 



^opt 



0.05 



/xl/10 



10" lO'^ 



10° 



Figure 3. (a) Final (steady state) average phonon number minimized with respect to 
f/m and [min{nj}] as a function of the thermal mechanical occupancy n{ujjn) [solid 
lines] . The lower bound rixL given by Eq. ( 54 1 is shown for comparison [dashed lines] . 
(b) Optimal normalized optomechanical coupling gj^'/w™ for which the minimum 
mechanical occupancies shown in (a) are attained. In both (a) and (b), Q„i = 5 x lO**, 
the light blue lines correspond to nj, = and n/ojm = 10, the dark blue lines to rip = 



and n/ujm = 0.1, 



and the red lines to K/ujm — 0.1 and a finite Up determined by 



T+'^™16 



with 



R 



(62) 



(63) 



It is straightforward to show that in the hmits g-m ^ with < 0, and <^ gm ^ 
with = —ijJm we recover, respectively, the approximate expressions (35), (53) [to 
lowest order in 7m/r]- 

The analytical formulas Eq. (61), (62) provide a quantitative basis for a precise 
analysis of the final temperatures that can be attained for arbitrary values of gmj In 
particular they can be readily minimized numerically with respect to the drive 's detuning 
and input power subject to the constraint imposed by the stability condition R~^ > 
that emerges from the Routh-Hurwitz criteria. This is the natural optimization problem 
that is posed by an experimental realization in which the natural linewidth of the cavity 
is fixed or hard to modify (admittedly sweeping this parameter within some range should 
be straightforward in electromechanical setups). We have performed this optimization 
using instead the exact expression that follows from Eq. (55), i.e. without the small 
7m approximation so that the transition to the trivial regime min{nf} = n{uJm.) is 
captured. Figure |3](a) compares this optimum (solid lines) as a function of the thermal 
equilibrium occupancy n{u!m) with the lower bound utl (dashed lines) furnished by 
Eq. (54) for representative values of the other parameters. The blue lines correspond to 
Up = (i.e. Up/ujm — > oo) for values of K/um in the Doppler regime (light blue) and in 
the RSB regime (dark blue). The red lines show the latter with Up/um = 10'^ instead. 
There are three distinct regimes: (i) n{ujm) < n^L so that the optomechanical coupling 
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Figure 4. (a) Final (steady state) average phonon number minimized with respect to 
gm and A^, [min{n^}] as a function of the thermal mechanical occupancy n{uj„i) and 
the ratio n/ajm, for Q„i — 5 x 10^ and n-p = 0. (b) Idem as a function of the thermal 
mechanical occupancy n{uj,n) and the ratio ajp/cum, for Qm — 5 x 10^ and n/LOm = 0.1. 



only raises the temperature and g'^^ = 0, (ii) the minimum occupancy is determined 
by the fundamental quantum backaction limit and is thus independent of the ambient 
temperature [i.e. constant as a function of n{uJm)], and (iii) the minimum occupancy 
is determined by classical noise and is thus linear in the ambient temperature. In this 
latter regime for finite Up/um there is a noticeable transition at which the minimum 
deviates from the Up = result as the cavity becomes thermally activated. One can also 
note that in this linear regime utl is only reachable in the RSB limit. In turn Fig. |3](b) 
shows the corresponding optimal values for g^- The saturation at high temperatures 
just reflects the fact that (n™ + 'nFj)/n{(jJm) becomes temperature independent in that 
limit and attains its minimum at a finite gm [cf. (62)]. Finally, Fig. |4]^a) and (b) 
illustrate, respectively, the dependence of the optimum with the ratios K/oJm and ujp/ojm 
as the ambient temperature is varied in the RSB regime. The fundamental quantum 
backaction limit results in a distinct shoulder. In Fig. |4](a) its lower edge traces a 
line with a slope consistent with the quadratic dependence on njuj^. It is clear from 
Fig. |4](b) that for ratios Up/uj^ ^ 10^ and n(ujm) ^ 10^ (parameters that are relevant 
for electromechanical setups [21]) thermal noise in the cavity needs to be taken into 
account. 



5. Conclusion 



In summary, we have analyzed the cavity-assisted backaction cooling of a mechanical 
resonator in the quantum regime by deriving effective master equations for the relevant 
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degrees of freedom both in the perturbative and strong couphng regimes. These provide 
a description of the coohng dynamics that allows to establish a simple lower bound for 
the final occupancy [Eq. (54)] that can only be reached in the RSB regime. This bound 
implies that ground state cooling is only possible when the cavity linewidth is much 
larger than the heating rate induced by the mechanical dissipation but much smaller 
than the mechanical oscillation frequency, and the equilibrium thermal occupancy of 
the cavity is well below unity. In addition we give an analytical expression for the final 
occupancy valid in all regimes of interest that allows for a straightforward optimization 
of the parameters of the drive. Finally we analyze the dependence of this optimum on 
the ambient temperature, the cavity linewidth, and the ratio of the cavity's frequency 
to the mechanical frequency. 
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